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Abstract 



C ■ Carbon nanotubes have properties depending on the arrangement of carbon atoms on the 

tube walls, called chirality. Also it has been tried to grow nanotubes of only one chirality for 
more than a decade it is still not possible today. A narrowing of the distribution of chiralities, 
however, which is a first step towards chirality control, has been observed for the growth of 
nanotubes on catalysts composed of nickel and iron. In this paper, we have calculated carbon- 
metal bond energies, adhesion energies and charge distributions of carbon nanotube caps on 
Ni, Fe and NiFe alloy clusters using density functional theory. A growth model using the 
calculated energies was able to reproduce the experimental data of the nanotube growth on the 
alloy catalysts. The electronic charge was found to be redistributed from the catalyst particles 



to the edges of the nanotube caps in dependence of the chiral angles of the caps increasing the 
reactivity of the edge atoms. Our study develops an explanation for the chirality enrichment in 
the carbon nanotube growth on alloy catalyst particles. 
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Introduction 



Carbon nanotubes have mechanical and electrical properties, depending on the helicity /atomic 
structure of the carbon atoms on their tubes surface (chirality). 1-3 The self assembling growth 
process of nanotubes generates nanotubes with various chiralities in the same growth ensemble. 4 
Three phases have been suggested to account for the chirality distribution found in nanotube en- 
sembles, the nucleation phase, 5 the growth phase, 6 ' 7 and the termination of the growth. 8 The 
nucleation phase involves the formation of graphene-like sheets (sheets of hexagonally arranged 
carbon atoms) which need to have pentagons in their structure to induce curvature. 9 Only with the 
curvature, which is achieved by a curved template particle, it is possible for a graphene sheet to 
transform into a carbon nanotube cap and lift off the particle.—^ This cap can then elongate with 
a chirality dependent growth rate by carbon addition to the edge of the growing tube.-^ 7 - 1 ^ 2 - Many 
theoretical and experimental studies tried to understand how the self assembly of nanotubes works, 
however, it is still not fully understood how the chirality can be controlled during the growth of 
carbon nanotubes. £&&12rm 

In recent experimental studies it was found that growth on alloy NiFe particles leads to an 
enrichment of certain chiralities.— ^ A theoretical study on a nickel particle suggested that the 
electronic charge transfer might be important to control the chirality-selective growth process.— 
Another theoretical study attempted to explain the chirality enrichment on the NiFe alloy particles 
by the lower excess energies for certain chiralities. 23 The study, however, considered plane surfaces 
and elongated nanotube caps (non-minimal seed caps with only an inferior number of growth 
sites/kinks), which does not seem to be appropriate to describe the chirality selection on a catalyst 
particle.— ) 2 4 - 

In this paper we study the (5,5) armchair and (9,0) zigzag carbon nanotube caps connected to 
various Ni, NiFe and Fe catalyst particles. We calculate the average carbon-metal bond energies, 
adhesion energies and charge distributions using density functional theory. We find the highest ad- 
hesion and lowest excess energy for armchair and zigzag caps on the Ni27Fe28 alloy cluster. Small 
energy differences between the armchair and zigzag bonds allow to derive carbon addition barriers. 
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The barriers can be used in a growth model leading to chirality distributions that compare satis- 
factorily with the recent experimental results of growth on NiFe alloy particles. A charge transfer 
from the catalyst particles to the caps induces a dipole moment between the catalyst particle and the 
cap. The polarity of the bond between the cap edge and catalyst atoms increases with increasing 
Fe content. The charge redistribution is found to depend on the chirality, as the line density of edge 
sites increases with higher chiral angle. The excess electron charges on the armchair rim atoms are 
found to be between (2.90 ± 0.06) e for Ni and (4.15 ± 0. 14) e for Fe. The increasing reactivity 
induced through the excess electron charges on the edge carbon cap atoms allows to explain why 
the nanotube growth rate on iron is higher than on nickel. 

Methodology 
Computational Methods 

We performed spin polarized density functional theory calculations with the ab-initio package 
SIESTA. 25,26 We used the generalized gradient approximation parameterized by Perdew, Burke 
and Ernzerhof, 27 as the bias towards compact cluster structures is reduced compared to the local 
density approximation. 28 The calculations used norm conserving nonlocal pseudopotentials gen- 
erated in the scheme of Troullier and Martins with the parameters presented in the PhD thesis 
of Dumlich. 24,29 The valence electrons were described by localized pseudoatomic orbitals. To 
balance the computational time and the accuracy to a reasonable level, we used a double- £ po- 
larized (DZP) basis set. The cutoff radii of the orbitals were chosen with r s = 6.099 Bohr and 
r p = 7.832 Bohr for the s and p orbital of the carbon atoms, r s = 9.649 Bohr and rj = 6.001 Bohr 
for the s and d orbital of the iron atoms, and r s = 9.187 Bohr and rj = 5.572 Bohr for the s and d 
orbital of the nickel atoms. The mesh-cutoff for the real-space integration corresponded to about 
350 Ry. We used only the T-point to calculate the total energies, as all studied systems have finite 
dimensions. 

For our calculations we consider the situation, in which a carbon nanotube cap has already 
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formed on a catalyst particle, but is not elongated. The systems studied in this paper therefore 
consist of two parts, a catalytic particle and a carbon nanotube cap. We consider the Fess, Nii2Fe43, 
Ni27Fe28, and N155 clusters to understand the influence of alloy systems and chiralities on the cap- 
cluster interaction. 




Figure 1: Top view on a ball and stick sketch model of the spots of a 55 atom (iron) cluster. The 
atoms that form the top of the cluster are marked, a) Spot 1, fcc(lOO) with a 3 times 3 atom square 
(9 atoms), b) Spot 2, fcc(lll) with a 3 atom sided triangle (6 atoms), c) Spot 3, with only 1 atom 
at the top of the cluster. 

Our (deformed) icosahedral catalyst particles consists of 55 atoms 30 that initially form a higly 
symmetric structure containing six fcc(lOO), eight fcc(l 1 1) surfaces, and three distinctive spots to 
add a carbon nanotube cap, see[T]a) to[T]c). The icosahedral particles have been found to be the 
most stable configuration. 30 The relaxed particle structures do not keep the icosahedral symmetry 
and are not expected to represent the global minimum, as it is not feasible to find the global min- 
imum of the clusters in density functional theory.— ^ The energy differences between different 
geometry optimized catalyst clusters range between 2 eV and 3 eV, corresponding to previously 
reported values.— We calculated the energies of the geometry optimized catalyst structures from 
each individual initial structure of the combined system to derive consistent carbon-metal and ad- 
hesion energies. 

The carbon nanotube caps were created with the program code CaGe using the isolated pen- 
tagon rule (IPR), which states that caps are energetically most stable, if all six pentagons needed for 
the cap inclination are isolated from each other.— - — We created and geometry optimized fullerene 
structures starting from the as-generated caps of CaGe to get to a decent cap structure. There are 
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many cap configurations for a certain chirality, which hinders to draw conclusions for the most 
chiralities.-22aM The isolated pentagon rule used for the generation of our caps, however, allows 
to reduce the number of possible caps, by only considering energetically favorable structures. 34 
Especially the (5,5) (30 atom cap) and (9,0) (39 atom cap) chiralities only have one possible cap 
structure which fulfills the isolated pentagon rule, 35 making them perfect candidates for our study. 



Figure 2: Ball and sticks sketch model of carbon nanotube caps on a Nii2Fe43 alloy cluster on 
three different spots. Carbon atoms in grey, nickel atoms in green, and iron atoms in orange, (a) 
Spot 1: (9,0) cap shows a Klein-edge through bond break of a pentagon at the rim. (b) Spot 2: 
(9,0) cap showing a hexagon and pentagon at the zigzag rim. (c) Spot 3: (5,5) cap showing an 
armchair rim formed exclusively by pentagons. 

The caps were transferred on the three different spots of the catalyst particles, see [TJ The fit 
of the rims of the carbon nanotube caps to the spots of the clusters were performed by hand. In 
[2] we show geometry optimized structures of the (9,0) and (5,5) cap bound to the three spots of 
the Nii2Fe43 catalyst as examples for all the systems we studied. The atoms at the rim of the 
cap bind to the cluster, see |2] The deformation (structure change) of the cluster is significant for 
the presented alloy systems. The atoms in the rim of the cap also adjusted their positions, which 
suggests a dynamic process of carbon and metal reshaping. 36,37 The edges/rims of the caps are 
composed of only armchair and zigzag sites, which is true for all nanotubes. If a pentagon bond 
gets broken at the edge, a Klein-edge can form and offer a site for addition of a single carbon atom 
to close the edge with a hexagon, see [2] a). The Klein-edge configuration occurs commonly for the 
(9,0) cap on the square spot and for the (5,5) cap on the top spot and might be a possible way 
to avoid the initiation barrier for a new layer. The configuration with the Klein edge was recently 
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suggested to be energetically favorable for the armchair growth of graphene (A5' site). 38 In general 
there is a nearly continuous number of ways to combine the cap with the cluster in dependence of 
the chirality of the cap, however, the chosen spots are expected to allow for the best comparison. 

All optimized geometries were relaxed to a maximal atomic force of 0.04 eV/A. The total 
energy of the combined system as well as the energy of the cluster and cap were calculated in 
various structures and basis sets to account for the basis set superposition error (counterpoise cor- 
rection).—"^ 2 The basis set superposition errors for our calculations range between 2.0 eV and 
3.0 eV with an average error of (2.5 ± 0.3) eV. The adhesion energy between the carbon nanotube 
cap and the catalyst cluster can be calculated with: 



■p z? -prelax prelax pBSSE /u 

^ad—^tot £<cap ^ cluster c w 



where E tot is the energy of the geometry optimized combined system of the cap and the catalyst 
cluster, is the energy of the geometry optimized cap, and E r c f™ er is the energy of the geometry 
optimized cluster.- 3 -^ 3 - 

To assess the stability of the combined system we compare the cap-on cluster system with other 
systems that contain the same number of carbon and metal atoms. This leads to formulas for the 
excess energy: 

K = Etot ~ E^ - ncE l c - E BSSE , (2) 

with E l c the energy per carbon atom for the system the combined cap-on-cluster system is compared 
to and iiq the number of carbon atoms in the cap of the combined system.-i 3 -^ We calculate 
the excess energy per atom in comparison to the energy of a system of an isolated metal cluster 
and a fullerene, to determine which system configuration is lower in energy and therefore more 
stable.-^ 3 -^^ Alternatively the excess energy compared to the fullerene structure can be obtained 
by removal of the dangling bond contributions of the adhesion energy. This leads to 

fullerene = ^ = ^ _ 2 . m . E ™ C - (n - m) • E™ , (3) 
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where E aa - is the adhesion energy, E™ c is the armchair bond energy in vacuum, and E™ c is the 
zigzag bond energy in vaccum. The factor 2 ■ m results from the number of armchair bonds at the 
rim of the nanotube cap and the factor (n — m) results from the number of zigzag bonds at the rim 
of the cap.- ' 45 The vacuum bond energies (for a straight cut rim) can be derived with the equations 

E c /2-E { c 5 a ' 5) E c /2-£^' 0) 

E^ ac = — 60 2-m C " P an( * Ez a ° = — 1% n-m ~ wnere ^c, is tne energy of a fullerene formed from 
two (n,m) caps. We used the system specific vacuum bond energies to determine the carbon-metal 
bond energies. To determine the carbon-metal bond energies per bond we divided through the 
number of bonds at the edge n + m. 

An error for the energies was estimated from the standard deviation by averaging over identical 
systems with small changes in their initial configuration. The standard deviations are rather large, 
as only two values were included which does not have a statistical significance, however, it allows 
to estimate the order of the error with at least 0.1 eV for the excess energy. The error for the 
adhesion energy is estimated to be slightly higher with about 0.7 eV. 

We performed Bader population analysis calculations on the cap-cluster systems to determine 
the electron charge transfer between the carbon cap atoms and the catalyst particle atoms.— 

Growth model 

The growth model applied in this paper is based on the idea that the carbon addition occurs to 
the hexagonal rim of the nanotube. 7,45 The hexagonal rim contains armchair and zigzag edges.- 14 - 5 - 
The addition of carbon dimers elongates the rim with hexagons and transforms the rim through a 
change of the number of energetically favorable addition sites. 7 ' 45 An energy barrier occurs for the 
addition of carbon atoms to the rim, which is especially high for the addition to zigzag rims.- 1 ^ 
The energy barrier for the rim transformation occuring through the carbon dimer addition to an 
armchair rim (initiation/closing of a new layer) is determined by 

A a = 2-\E z -E a \, (4) 
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where E a is the energy of an armchair carbon-metal bond and E z is the energy of a zigzag carbon- 
metal bond. 7 ' 24 The energy barriers influence the chirality dependent growth rate.-^ 7 - 1 ^ 2 - A factor for 
the different growth rates in dependence of the chirality is given with 



Ka.aa (n,m) ■ 8 a + A fla ., (n,m) ■ S a: _ ^ q 

r(,um)={ n+m ~ (5) 

A aa . z (n,m)-8 az otherwise, 



where 8 a = exp (—A a /kbT) is an exponential factor to account for the temperature dependence of 
the addition barrier to armchair sites and 8 az = exp (—A az /kt,T) = 1, as an addition barrier for aa.z 
sites (kinks) is negligible (A az = 0).^^ The other terms are defined as A aa aa = 2m — n — 1 + 
1 / (2m — n) and A aa z = min (m, n — m). 1 The growth rate factor T(n,m) can have values between 
and 0.5 in dependence of the chirality. 

r(n,m) allows to determine the relative difference in growth rate between different chiralities. 
To derive a chirality distribution we include the influence of the nucleation phase, i.e., whether a 
particular tube cap is nucleated or not. 34 We assume the tube diameters (and also the chirality) to 
be fixed by the nucleation. 5 ' 50 To consider the dependence of nanotube diameters on the diameter 
of the catalyst particles,— - — we therefore multiply Equation (??) by a Gaussian distribution of the 
nanotube diameters / (d; jU, a 2 ) and obtain a growth rate factor T* (n,m) which can be compared 
to a chirality distribution, 

1 (rf-M) 2 

r (n,m) = — i= e 2 ° 2 ' r («, m) . (6) 



The tube diameter distribution might also include additional effects that do not result from the 
particle diameters and which might not be covered by the Gaussian distribution, however, the 
distribution serves the simplicity of the model. 
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Results 



Adhesion energies and carbon-metal bond energies of nanotube caps 

The adhesion energies are high with slightly more than —20 eV, varying about a few eV in de- 
pendence of the spot, catalyst and cap. Therefore the nanotube caps are not expected to lift off 
the catalyst particle spontaneously and perform dome closure.-^^ The armchair caps have adhe- 
sion energies that are a few eV higher than the zigzag energies, which is a result of the number of 
dangling bonds at the edge of the caps. The zigzag caps have n + m = 9 + = 9 dangling bonds. 
The armchair caps have n + m = 5 + 5 = 10 dangling bonds. Therefore the adhesion energy of 
the armchair cap is higher than the adhesion energy of the zigzag cap, even though the energy per 
dangling bond is lower for the armchair compared to the zigzag bond.— 

Table 1 : Comparison of adhesion energies and carbon-metal bond energies of nanotube caps on 
metallic/alloy clusters in eV. The errors are standard deviations from averaging over the spots. 



cap 


particle 




E ad (eV) 


£ CM /bond (eV) 


Ecm (eV) 


(5,5) 


Fe 5 5 


(- 


-21.2±0.0) 


(0.43 ±0.00) 


4.3 


(5,5) 


Nii 2 Fe 43 


(- 


-23.2±1.1) 


(0.30 ±0.11) 


3.0 


(5,5) 


Ni 27 Fe 28 


(- 


-23.5 ±1.0) 


(0.27 ±0.07) 


2.7 


(5,5) 


Ni 55 


(- 


-23.2±0.3) 


(0.32 ±0.03) 


3.2 


(9,0) 


Fe 55 


(■ 


-20.7 ±-.-) 


(0.40 ±-.-) 


3.6 


(9,0) 


Nii 2 Fe 43 


(- 


-21.2±0.9) 


(0.33 ±0.10) 


3.0 


(9,0) 


Ni 2 7Fe 2 8 


(- 


-21.5 ± 1.1) 


(0.30±0.13) 


2.7 


(9,0) 


Ni 55 


(- 


-20.2 ±0.0) 


(0.46 ±0.00) 


4.1 



Averaging the adhesion energies for the alloy compositions we find Ejj 55 = (—21.0 ± 0.3) eV, 
£Ni 12 Fe 43 = (_ 22 3 ± L4 j eVj £ Ni 27 Fe 2 8 = (_ 223 ± L5 ) eV; and £ Ni 55 = ^_ 22A &y Two 

effects account for the adhesion energy. One is the structure and the other is the material. The 
mixture of two materials distorts the catalyst structure as they have different electronic structure 
resulting in different bond lengths and lattice constants. We observe an increase of adhesion energy 
for the NiFe alloy systems compared to the pure Fe clusters, see [TJ the effect, however, is not 
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significant. 

The adhesion energies averaged on the three spots introduced in Q3 show no obvious trend 
either, with E ad = (—22.0 ± 1 .5) eV for the square spot I, E^ d = (—23.0 ± 1 .4) eV for the triangle 
spot 2, and E\ = (—21.7 ± 1 .0) eV for the top spot 3. We use this result and average over different 
spots when calculating adhesion energies as a function of the cap and alloy composition, seeQ] The 
triangle spot shows the highest adhesion energy, but the adhesion energies calculated for the other 
spots are within the standard deviation of spot 2. 

The excess/carbon-metal bond energies decrease non monotonically with increasing Ni content 
from Fe to Ni. The Ni cluster shows the highest excess energy for the (9, 0) cap. The lowest excess 
energy is observed for the Ni27Fe28, which also showed the highest adhesion energy, suggesting 
that the Ni27Fe28 alloy cluster yields the best growth conditions for the systems compared in this 
study. The low excess energy allows for a fast formation of caps on the Ni27Fe28 cluster and the 
higher adhesion energy prevents the cap lift off, which makes the NiFe alloy systems more stable 
compared to the pure elemental catalyst clusters. 

Charge distribution on cap and cluster atoms 

Table 2: Charge redistributions between the carbon nanotube caps and the metallic/alloy clusters. 

is the total charge shift of the carbon atoms at the rim of the nanotube cap, qc is the total 
charge shift considering all carbon atoms of the cap, q? e is the total charge shift considering all 
iron atoms, and <?Ni is the total charge shift considering all nickel atoms of the catalyst particle. 



cap 


particle 


q™ (e) 


qc (e) 




<?Fe (e) 


qm (e) 


(5,5) 


Fe 55 


(4.15±0.14) 


(5.31 ±0.72) 


(- 


-5.32±0.72) 




(5,5) 


Nii 2 Fe 43 


(4.08 ±0.31) 


(5.06 ±0.52) 


(- 


-7.18±0.65) 


(2.12±0.46) 


(5,5) 


Ni 2 7Fe 2 8 


(3.67 ±0.24) 


(4.36 ±0.36) 


(- 


-6.46 ±0.20) 


(2.10±0.32) 


(5,5) 


Ni 55 


(2.90 ±0.06) 


(3.32±0.36) 






(-3.33±0.38) 


(9,0) 


Fe 5 5 


(3.26±-.-) 


(4.68 ±-.-) 


( 


-4.68 ±-.-) 




(9,0) 


Nii 2 Fe 4 3 


(3.39±0.17) 


(4.62 ±0.22) 


(- 


-6.71 ±0.42) 


(2.09 ±0.59) 


(9,0) 


Ni 2 7Fe 2 8 


(2.99 ±0.05) 


(3.96 ±0.02) 


(- 


-6.48 ±0.29) 


(2.52 ±0.28) 


(9,0) 


Ni 55 


(2.34 ±0.08) 


(3.11 ±0.42) 






(-3.11 ±0.42) 
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The charge population on the atoms in the rim of the cap and for the metal atoms in the catalyst 
particle shows a stronger trend than the energies, see |2] All charge values are excess charges 
compared to the isolated atom valence electron configuration with 8 charges on each Fe atom, 
10 charges on each Ni atom and 4 charges on each C atom. The metal atoms partly loose their 
electrons to the carbon atoms in the cap, with whom they form carbon-metal bonds. The amount 
of electron charge transfer to the carbon atoms depends on the catalyst element and on the bond 
type of the edge atom (zigzag or armchair). The charge on the rim of the armchair cap is higher 
than on the rim of the zigzag cap. The trend weakens, but does not vanish, if the charge transfer 
per bond is considered, as the number of carbon-metal bonds is 10 for the armchair and only 9 
for the zigzag cap. Considering the charge transfer to the whole cap and dividing through the 
number of bonds leads to an equal charge transfer to the armchair and zigzag caps per bond. The 
highest layer of the catalyst atoms supplies the major part of the electrons to the carbon nanotube 
cap. The carbon atoms of the cap that are not part of the rim have an average valence charge of 
about 4 e, with low deviations (below 0.1 e), which means that they do not take part in the charge 
redistribution/polar binding between the metal catalyst and the carbon cap. The electron charge 
redistribution is localised at the outer rim atoms of the cap that form the carbon-metal bonds, 
corresponding to a polar bond. The localisation is slightly higher for the armchair edges with 
about 80% of the charge localised at the outer rim atoms compared to about 74% at the outer edge 
zigzag cap atoms, which is the reason why the armchair rim atoms yield a higher charge per bond. 

The charge on the cap increases with Fe content and becomes maximal for the elemental Fe 
cluster. We observe the same behaviour at the rim, with the exclusion of the (9,0) cap on Fe, 
which, however, may be an artifact, as only one cap was considered on spot 3. The higher charge 
points to a higher reactivity on Fe compared to Ni, which likely leads to a faster growth rate on 
iron compared to nickel, previously observed experimentally. 56 

The Fe atoms do not only supply their electron charge to the C atoms, but also to the Ni atoms, 
see the charge transfer for the NiFe alloy systems in |2] The iron atoms loose about 5-7 electron 
charges. The nickel atoms either gain about 2-2.5 e in the alloy systems or loose about 3 e if Ni 
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is the only element in the catalyst particle. The carbon atoms in the cap always receive electron 
charges, see|2] The total amount of charge received by the carbon cap is about 3-5 e. The amount 
of charge supplied to the rim is slightly lower with about 2-4 e. 

Following from the reactivity argument, which results from the increased charge on the edge 
atoms, we can give a geometric argument for the preference of armchair over zigzag structures in 
the following. It follows from the number of edge sites. We consider the line density of edge sites, 
which corresponds to the number of edges sites divided by the circumference of the nanotube rim 



X = 



N a + N z n + m 



C h 



ao ■ Vn 2 + nm + m 2 



(7) 



Ch 



the circumference of the tube. The 



with N a + N Z the number of armchair and zigzag sites and 
righternmost equation follows for straight rim configurations. Equation (??) leads to a line density 
of X a = 2/ (\/3ao) for (n = m) armchair and X z = l/ao for (n integer, m = 0) zigzag tubes, for all 
other tubes (« / m ^ 0) the value of the line density is between X a and X z (X a > X c > X z ). Con- 
sidering the fact that a higher density of edge sites increases the number of carbon metal bonds, 
directly gives the argument why armchair tubes are prefered compared to zigzag tubes, as the num- 
ber of electrons at the edge is increased, yielding a higher reactivity. The comparison of the (5, 5) 
cap to the (9,0) cap presented here, serves as an example as they have nearly identical diame- 
ters, but a different number of edge sites. The increased number of sites with higher chiral angles 
leads to an enrichment of armchair/near- armchair tubes in nanotube samples from simple geomet- 
ric reasoning. This geometric argument might also be translated to growth rate considerations of 
nanotubes. The edges of tubes with higher chiral angles contain more sites for carbon atoms to 
dock, independent of the exact addition mechanism, supporting arguments for faster growth of 
higher chiral angle tubes. 

The charge transfer also induces a dipole moment in the nanotube cap and cluster system, see|3] 
The highest electric dipole moment can be found with 14.7 Debye for the (9, 0) cap on spot 3 of the 
Fe cluster and the lowest electric dipole moment is 7.4 Debye for the (5,5) cap on the Ni cluster. 
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Table 3: Electric dipole moments A between the nanotube caps and metallic clusters in Debye. 
The last column of the table shows electric dipole moments for the armchair/zigzag caps averaged 
over the spots of a specific catalyst composition. 



cap 


particle 


A avg (Debye) 


(5,5) 


Fe 55 


12.4 


(9,0) 


Fe 5 5 


14.7 


(5,5) 


Nii 2 Fe 43 


10.6 


(9,0) 


Nii 2 Fe 43 


12.1 


(5,5) 


Ni 2 7Fe 2 8 


10.0 


(9,0) 


Ni 27 Fe 28 


11.4 


(5,5) 


Ni 55 


8.6 


(9,0) 


Ni 55 


9.7 



A decrease of the electric dipole moment is correlated with the Ni content in the composition 
of the catalyst particle, where higher Ni content leads to lower electric dipole moments. The 
zigzag caps have higher electric dipole moments compared to the armchair caps, which is likely 
a result of the weaker localisation of electron charge at the carbon-metal bond forming atoms 
for the zigzag caps. The difference of the average electric dipole moments of a certain catalyst 
composition between armchair and zigzag caps decreases with increasing Ni content, with the 
highest difference of 2.3 Debye on Fe and the lowest difference of 1.1 Debye on Ni. The electric 
dipole moments generate an electric field, which was suggested to increase the landing probability 
of carbon atoms on the catalyst particle. 57 Therefore the landing probability on Fe is higher than 
on Ni, leading to a higher growth rate on catalysts containing Fe. A higher landing probability 
might also become a problem, if the carbon atoms cannot be incorporated at the edge faster than 
the carbon precursors land on the particle, as this might lead to amorphous carbon encapsulation 
of the catalyst and prevent further growth. The higher landing probability for zigzag caps can lead 
to a reduced number of zigzag tubes, as zigzag edges are expected to have a slower growth rate, 
which makes it especially hard to fulfill the requirement of fast carbon incorporation to prevent 
catalyst encapsulation. 
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Discussion 



In the following we want to discuss the results we can derive from the energies and charge dis- 
tributions calculated in this paper and compare them to other studies. There are two important 
results, that can be derived from the carbon-metal bond energies. First, the energies are important 
for the nucleation phase, as lower excess (carbon-metal bond) energies point to a higher forma- 
tion probability.-^^ Second, the carbon addition barrier A a can be derived from the carbon-metal 
bond energies for the growth phase, which was found to influence the chirality dependent growth 
rate. 6,7,12 Both phases determine the chirality of the carbon nanotube ensemble which is grown 
during the nanotube synthesis. We determine the barrier energies A a (Equation (??)), by using 
E a = Eq^ /bond and E z = E^^ /bond of Q] The barriers result from small deviations between 
armchair and zigzag bond energies. Previous studies considered flat metallic surfaces as catalysts 
and tried to optimize the fit between the catalyst surface and the edge of the cap. 5,23 On a curved 
particle, as in this paper, a perfect fit between the edge of the nanotube cap and the catalyst particle 
is not possible, which increases the carbon-metal bond energies. The bond energies for armchair 
edges from Reich et ah, derived for a flat Ni surface, range from E a = 0.12 eV to E a = 1.12 eV, 
comparing well to our average value E^ 155 = (0.32 ± 0.04) eV.— The values for zigzag edges from 
Reich et al. range from E z = 0.16 eV to E z = 1.44 eV, which are also comparable to our value 
E^ 155 = 0.46 eV. 5 The caps connect to various spots on the catalyst clusters which increases the 
deviation of the bond energies and therefore renders the chiral selectivity by structural fit of the 
cap and the catalyst particle unlikely. The large standard deviations point to a general problem 
for ab-initio studies of carbon nanotube growth. The quantitative reproducibility is rather weak 
and it should be desired to test various systems with slightly different configurations/parameters. 
The energy barriers for the studied catalyst compositions are calculated with A" 55 = 0.06 eV, 

A Nii 2 Fe 43 = q Q6 e y^ A Ni 27 Fe 28 = q and A Ni 55 = q 2g e y baniers ar£ equal for aU irQn 

containing catalyst compositions. 

The barriers can be inserted in Equation (??) to determine the chirality distributions, where 
the contribution of T (n, m) (Equation (??)) is equal for all systems that mainly contain iron atoms 
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Figure 3: Normalized abundances in dependence of the tube diameter; a) (jU = 9.0 A, o = 0.6 A) 
and c) (/! = 8. 1 A, o = 0.6 A) estimated from our theoretical growth model, see Equation (??); b) 
and d) from experimental photoluminescence data by Chiang et al. for nanotubes grown on Ni and 
on a nickel-iron alloy at 600° C— 
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(8 a = 0.45 for 600° C) and only the gaussian diameter distribution factor changes the chirality 
distribution. We can compare the results of our model to the experimental results presented by 
Chiang et ah, who grew nanotubes on NiFe alloy systems to analyse the influence of the catalyst 
composition on the chirality distribution of a nanotube ensemble.— ^ They derived chirality dis- 
tributions from photoluminescence data using calculated photoluminescence intensities.- 2 ^^ The 
chirality distributions satisfactorily fit the results derived from our growth model, see |3] 20 Espe- 
cially the growth on iron containing systems is interesting, as the experimental study found the 
chirality distributions grown from catalyst with the composition Nio.27Feo.73 and Nio.sFeo.s to be 
almost identical, 20 corresponding well to our results. The Ni catalyst particle has a significantly 
larger barrier energy, leading to a suppression of armchair growth sites, as, e.g. 8 a = 0.024 for 
600° C. In the experimental study the nanotubes grown on a Ni catalyst show a relatively wide 
chirality distribution with a peak for (9,4).— ^ Especially important seems to be the diameter 
region of the nanotubes with 9.0 A for (9,4) and 8.8 A for (7,6) which have the highest inten- 
sity/abundance in the experimental study. 20 A slight descent of intensity occurs for chiralities with 
diameters that have smaller/higher tube diameters than about 9.0 A, pointing to a lower number 
of catalyst particles, or other unknown effects, to grow tubes of that higher/lower diameters. In- 
creasing the iron content of the composition of the catalyst particles until iron becomes the major 
component, leads to a significant narrowing of the chirality distribution to only a few chiralities 
at lower diameters, compared to the Ni catalyst particle, see[3]d). The fcc-lattice constant of iron 
ap e = 3.45 A is lower than the lattice constant of nickel a^i = 3.63 A, 60 which leads to smaller 
diameter alloy catalyst particles with increasing Fe content, which might be a reason for the shift of 
the chirality distribution to lower diameter nanotubes ((7,6) and (8,4) abundance increased), as the 
diameter of the grown nanotubes depend on the diameter of the catalyst on which they are grown 
in the tangential growth mode/under growth conditions close to thermodynamic equilibrium. 
The experimental study, however, tried to obtain equal particle diameters through the preparation 
process, pointing to a dependence on the material instead of the catalyst diameters.— Our model 
was not intended to perfectly reproduce all abundances, however, it still gives a fair approximation 
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to the experimental results and successfully reproduces the most significant change in the chirality 
distribution by the change of the catalyst. 

Besides the chirality distributions obtained from our growth model we also want to compare 
our charge transfer results to the literature. Wang et al. suggested that the short ranged charge dis- 
tribution on nanotube edge atoms and catalyst atoms might be important for the chirality- selective 
growth of carbon nanotubes, as electron charges would increase the reactivity of the edge atoms.— 
We observe an increase of electron charge on the carbon edge atoms with charge supply by the 
metal atoms in agreement with Wang et al.. 22 The average charge values on the carbon rim atoms, 
see [21 compare well to the values calculated by Wang et al. for, e.g., the (5,5) nanotube cap on 
nickel we find an average value of 0.29 e which compares to the slightly higher values of Wang et 
al. between 0.31 e and 0.38 e. To determine the effect of the charges on chirality distributions, 
we also put a focus on alloy systems. In alloy systems a charge distribution between two metallic 
species leads to an electron accumulation not only on the carbon edge atoms, but also on the nickel 
atoms. We find higher charges on armchair than on zigzag edges, which was suggested by Wang et 
al. to be used to influence the chirality. 22 The alloy composition has a significant effect on the 
charge distribution. We find an increase of electron charge on the carbon cap edge atoms from 
Ni to Fe with increasing Fe content in the alloy, pointing to an increased growth rate of nanotubes 
through increased reactivity of the nanotube edge atoms, which compares well to the higher growth 
rates found for iron compared to nickel.— Another relevant factor for the growth rate was found to 
be the metal d orbital energy. 61 The charge distribution patterns suggested by Wang et al. resemble 
the edge structure of armchair and zigzag sites. 

Theoretical studies can only model some aspects of the nanotube growth, neglecting other as- 
pects, e.g., the effect of Ostwald ripening, 8 that influence the chirality distribution as well. Further 
the chirality distributions determined in experiments have to be regarded with care as huge differ- 
ences for the abundances of the chiralities might arise through the method used to determine the 
abundances, i.e. the intensity of a measured entity is not directly proportional to the abundance 
of the tube. 5 ^^ Therefore we did not expect to obtain results that perfectly match our growth 
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model, however, we see it as a success that the model correctly describes the qualitative features 
of the chirality selective growth process, which suggests, that the model might include some part 
of the truth to solve the puzzle of chirality selective growth. 

Conclusions 

In summary we calculated adhesion energies, excess energies, and electronic charge redistributions 
between carbon nanotube caps and NiFe alloy systems using density functional theory. The highest 
adhesion energies and lowest excess energies are found for the Ni 2 7Fe 2 8 alloy cluster, for both 
armchair and zigzag caps. The energy differences between armchair and zigzag were found to be 
low. The curved form of the catalyst particle can be regarded as a constraint to the fit between 
the nanotube edge and the catalyst, which tends to lower the energy difference between armchair 
and zigzag caps. The small energy differences between the armchair and zigzag caps allow to 
derive carbon addition barriers, which - using the growth model presented in this paper - lead to 
chirality distributions that compare satisfactorily with experimental results. The charge transfer 
between the cap and the catalyst particles increases with increasing Fe content, which induces a 
dipole moment. The charge transfer to the armchair caps is higher than to the zigzag caps, in 
contrast to the electric dipole moment, which is higher for zigzag than for armchair caps and has a 
maximum of about 15 Debye on the iron particle. This is a consequence of the chirality dependent 
line density of edge sites, which decreases with lower chiral angle. The excess electron charges 
on the carbon rim atoms increase with Fe content of the catalyst particle from (2.90 ±0.06) e for 
Ni to (4. 15 ± 0. 14) e for Fe. The excess electron charges increase the reactivity of the carbon cap 
atoms, which explains why the nanotube growth rate on iron is higher than on nickel. Our results 
will be useful for the understanding of the growth of carbon nanotubes on alloy catalysts. 
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